Screening in YBa 2 Cu 3 07_5 at large wave vectors 
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We present experimental inelastic x-ray scattering (IXS) and ab initio time-dependent density- 
functional-theory (TDDFT) studies of YBa2Cua07_5. The response of the low-lying Ba 5p and Y 
4p core electrons is shown to interact strongly with the Cu 3d and O 2p excitations, with important 
consequences on screening. The agreement between IXS and TDDFT results is excellent, apart from 
a new type of excitations, mainly related to loosely bound Ba electrons and significantly affected 
by correlations. This points to correlation mechanisms not fully described by TDDFT that might 
have a role in giving rise to antiscreening. 



I. INTRODUCTION 

Despite 20 years of research since the discovery of 
cupratcs as high-temperature superconductors (HTSC)fi 
the nature of the superconductivity and the pairing 
mechanism in these systems still remain unknown. Ex- 
tensive studies of HTSC cuprates have been done using 
the most varied techniques and on almost every observ- 
able imaginable. The dielectric function e(r, r',t — t') is 
in particular a key quantity for superconductivity. Its in- 
verse e _1 measures the screening of the bare Coulomb re- 
pulsion and can directly indicate reverse screening (anti- 
screening) spatial regions where the interaction between 
two electrons is attractive rather than repulsive. In these 
regions electrons pair up, giving rise to superconductiv- 
ity. The study of can thus unambiguously indicate 
such domains and shed light to the pairing mechanism. 

The dielectric function is probed often via its Fourier 
transform e(q, ui) as in reflectance and cllipsometry spec- 
troscopy, but usually only at negligible transferred mo- 
mentum q| = q — >• 0. This provides access only to long- 
range screening (|r — r' — > oo). In HTSC the coher- 
ence length is of the order of few inter-atomic distances 
[e.g. 1.5 nm in YBa 2 Cu 3 07_<5 (YBCO) in the ab plane^]. 
The antiscreening domains have to exist at comparable or 
shorter distances. In the large- q range, which would cor- 
respond to short-range correlations, experimental tech- 
niques are scarce due to the high-kinctic-energy particles 
required as a probe. As a consequence, this dynamic do- 
main of £ _1 (q, w) for q > has been a no-man's land 
with very few applicable experimental probes. There 
are basically only two techniques available: the electron 
energy- loss (EELS) (Refs. i and i) and inelastic x-ray 
scattering (IXS) (Ref. H) spectroscopies. Both, in prin- 
ciple, measure the loss function, which is proportional 
to — Im [e~ 1 (q, w)] . EELS, however, is inherently lim- 
ited to relatively small values of q due to an increasing 
multiple scattering at large wave vectors. IXS is left as 
the only technique probing exchanged wave vectors that 
correspond to inter-atomic distances, the typical range 
accessible to IXS being q > 5 nm -1 . It can hence be a 



valuable technique to study HTSC. 

IXS studies in strongly-correlated materials have been 
so far limited by the low count rates in high-Z elements. 
Thanks to the rapid development in instrumentation, it 
is now becoming possible to apply IXS even to HTSC 
cuprates. Resonant IXS (Refs. has already been 

used to study such systems but these experiments do not 
probe the dielectric function and should not be confused 
with non-resonant IXS which we use here. In particular 
on YBCO, EELS experiments^— were performed soon 
after the discovery of its superconductivity, but there are 
no existing IXS experiments. This work is an attempt to 
close the gap in our accessible kinematic probing range 
of excitation energies and momenta. 

The dominant theoretical picture on HTSC relies on 
the Hubbard modeli 3 - and invokes a strong correlation 
mechanism to explain the pairing and the superconduc- 
tivity in cupratesi 1 ^ This model accounts only for the last 
unpaired electron of a copper atom and disregards all the 
complicated chemical, atomic and electronic structures of 
HTSC cuprates. They are supposed only to renormalize 
the two adjustable parameters of the model. Today, a 
good framework to search for a solution to the Hubbard 
model is the dynamical mean-field theory (DMFT) (Ref. 
[THh and more recently cellular DMFT which introduces a 
k dependence necessary to address the evident anisotropy 
in HTSCi^Li 7 . However, systematic improvements within 
this approach are difficult without introducing new ad- 
justable parameters. 

It has been recently demonstrated that ab initio theo- 
ries can be applied to HTSC. Density-functional theory 
(DFT) has show n 18 i 19 its ability to reproduce the ground- 
state atomic structure of YBCO within its typical error 
< 3%. However, DFT does not describe excited states. 
On the other hand, time-dependent density-functional 
theory (TDDFT) (Refs. HO and HO) is in principle an 
exact theory to describe neutral excitations and the di- 
electric function. In practice, since the exact exchange- 
correlation functional of TDDFT is not available at the 
moment, approximations such as the adiabatic local- 
density approximation (TDLDA) (Ref. l22t ) are required. 
The latter has shown to provide loss-function spectra 
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FIG. 1: Experimental and calculated energy-loss functions of YBCO for (a) q || c and (b) q _L c. The most prominent spectral 
features are labeled A-G and highlighted with dashed lines in the experimental plots. The nominal core-electron excitation 
energies are marked with dashed lines and labeled in the theory plots. For clarity, the spectra are shifted vertically from each 
other by an amount proportional to the corresponding momentum transfer value. Note that the experimental S(q,Lj) are 
area-normalized to the TDLDA theory as described in the text. The experimental geometries for the two directions of q are 
shown on the left. 



for ordinary semiconductors and insulators, e.g. silicon, 
in good agreement with both EELS (Ref. (H) and IXS 
(Ref. l24l ) experiments. Recent works have shown that 
TDLDA provides results in excellent agreement also with 
the exact solution of a one-dimensional Hubbard model, 
reproducing spin and charge collective excitations in a 
Luttingcr liquid i 25 ' 26 TDLDA can hence be a good for- 
malism even in strongly correlated systems. However, so 
far there are no TDDFT studies of the dielectric function 
in cuprates. 

In this work we present an investigation of the in- 
verse dielectric function e _1 (q, w) on a prototypical high- 
temperature superconductor cuprate, YBa2Cu307_<5 
(6 = 0.07) combining IXS experiments with TDDFT 
calculations. IXS probes the dynamic structure factor 
5(q, u), which is related to the dielectric function via 

S(q,w) = -(% 2 )/(47r 2 e 2 n) Im £ _1 (q,a;) (1) 

n being the electron density. We study S(q || c,ui) and 
5(q _L c,w) for energies of 5-60 eV and q between 5.6- 
30.0 nm^ 1 . The comparison of experimental results with 



TDDFT calculations provides an important benchmark 
for studying the electronic response of cuprates in the 
whole physically relevant momentum-transfer range. In 
particular, the large energy range studied here together 
with calculations facilitate the identification of the role 
of yttrium and barium in the valence electron dynamics 
of YBCO. 



II. EXPERIMENT 

The non-resonant IXS measurements of the loss func- 
tion of YBCO were done at the beamline ID16 (Refs. H?] 
and |28| ) of the European Synchrotron Radiation Facil- 
ity (ESRF). The measurements were performed at room 
temperature using monochromatic x rays with energies 
of 7.9-8.0 keV. The spectrometer was based on a spheri- 
cally bent Si(444) analyzer crystal in the Rowland-circle 
geometry. The bending radius of the crystal was 1 m 
and the Bragg angle was fixed to 89° . The spectrometer 
observed the intensity of scattered photons with a fixed 
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energy, and the energy transfer was tuned by changing 
the incident-photon energy. The sample was in a shape 
of a plaquctte with a size of 5x5x0.1 mm 3 , with the 
polished main face oriented perpendicular to the c axis 
and the edges along the a and b axes. The sample ori- 
entation was analyzed using a Laue photograph. The 
superconducting transition of the sample takes place at 
Tc = 94 K. The measurements were performed with q 1 1 c 
and with qlc. The measurement geometries in the two 
cases are shown in Fig. [T] While the sample is easy to 
align unambiguously for the measurements with q || c, 
the measurements with qlc were slightly more dim- 
cult with the sample in question. First of all, the usual 
twinning of a YBCO was present in the sample, giving 
a small amount of uncertainty to the exact directions of 
the a and b axes. While the measurements were per- 
formed as close to q | a as possible, for this reason we 
refer to these results rather as q _L c. Furthermore, the 
measurement in the ah plane was done in a near-grazing 
incidence and near-grazing exit geometry. The incidence 
and exit angles were thus small but finite, giving a small 
c component to the momentum-transfer vector as shown 
in Fig. [IJb). However, this component was kept always 
smaller than the q resolution and thus is expected to have 
a negligible effect on the results. 



All spectra at each fixed q were measured several times 
in the same conditions and after a normalization to the 
incident-beam intensity were found to be identical within 
statistical accuracy. The resulting spectra were then av- 
eraged in order to increase the statistical accuracy. This 
procedure minimizes the possibility for any experimen- 
tal instabilities during the measurements. The resulting 
spectra, which at this stage are scattering intensities as 
a function of energy transfer uj and momentum trans- 
fer q were corrected for sample self- absorption and linear 
background. We refer to the Q = as the zone center, 
and thus momentum transfer values given below are ab- 
solute. The energy resolution of the experiment was 
1.0 eV and the momentum-transfer resolution 1.5 nm _1 . 
The zero-energy-loss quasielastic line was subtracted by 
using an exponential fit to its positive-energy-transfer tail 
above 2.5 eV. This procedure yields reliable spectra above 
~ 3 eV. The experimental spectra were normalized to the 
same area with the theoretical TDLDA spectra. Since 
in the current computational approach 2 ^ the theoretical 
results are limited to values q = where n is an in- 
teger and x the length of the direct-space lattice vector, 
it was not necessarily possible to perform the calcula- 
tion of S"(q, to) exactly at all measured values of q. For 
the determination of the area-normalization factor, we 
first determined numerically the theoretical behavior of 
the integrated area S^ OTy (q) = J 6 ° cV S(q, lj)Auj and in- 
terpolated the resulting curve to the values of q which 
were actually measured. A normalization to this area- 
integrated value yields finally the experimental S*(q, o>). 



III. THEORY 

Our computational starting point was a standard 
ground-state DFT LDA calculation (as in Ref. [HI) of 
the total energy and the electronic density in YBCO 
using the ABINIT code.— A plane-wave basis set was 
used with a converged 160 Ry kinetic-energy cutoff 
and periodic boundary conditions. As turns out later 
to be very important, we included the low- lying-core 
states Ba 5s5p, Y 4s4p, Cu 3s3p, and O 2s in the va- 
lence using Hartwigsen-Goedecker-Hutter pseudopoten- 
tials. Employing the calculated DFT Kohn-Sham ener- 
gies and wave functions we carried out a TDDFT cal- 
culation of the polarizability and the dielectric function 
using the DP code^ using both TDLDA and random 
phase approximation (RPA). The latter neglects corre- 
lation effects in the dielectric response of the system, 
and the comparison of TDLDA and RPA might give a 
hint where correlation effects play an important role. 
The independent-particle polarizability was constructed 
summing over 350 Kohn-Sham bands calculated at a 
Monkhorst-Pack 8x8x3 fc-point grid shifted by \ 
The dimensions of polarizabilitics and longitudinal 
dielectric functions as matrices in the reciprocal space 
were 347x347. The inversion of the polarizability matrix 
produces so-called local field effects where different ex- 
citation channels mix. This effect, usually neglected in 
the localized description, is well known from earlier first- 
principles studies on, e.g., transition metals.— It should 
be emphasized that we found these local-field effects to 
be crucial to obtain a good agreement between experi- 
ment and theory and should definitely not be neglected 
in this system. We used experimental lattice parame- 
ters, and we verified that the relaxed theoretical atomic 
structure^ produces only minor changes to the spectra. 
All theoretical data shown here were convoluted with the 
experimental resolution function. 

IV. RESULTS AND DISCUSSION 

The measured and calculated spectra are shown in Fig. 
[1] The agreement between the experimental and theoret- 
ical results is remarkably good. This is a manifestation 
that the TDDFT is an extremely useful tool in under- 
standing the physics of even this canonical strongly corre- 
lated system. In fact, all relevant physics of the electronic 
excitations studied here are captured by the TDDFT. 
The energy range of the investigated excitations covers 
the collective oscillations of the valence electrons (plas- 
mons), consisting mostly of the hybridized Cu 3d and O 
2p, and the excitations of shallow core electrons of Y and 
Ba. Because the latter interact strongly with the valence 
electrons due to the important overlap in energy, we call 
them semicore states. The literature value d 32 ' 33 for their 
atomic excitation energies are 14.8 eV (Ba 5p3/ 2 ), 17.0 
eV (Ba 5p 1/2 ), 23.1 eV (Y 4p 3/2 ), 24.4 eV (Y 4p 1/2 ), 30.3 
eV (Ba 5s), 41.6 eV (O 2s) and 43.8 eV (Y 4s). These are 
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FIG. 2: Zoom- up to the region of 20-32 eV energy transfer to 
show the Ba 5p + Y 4p related dispersing excitations C-E in 
experiment, TDLDA and RPA. Their dispersion are marked 
with dashed lines as a guide for eyes. 

marked in the Fig. Q] with dashed lines. Assuming that 
these electrons indeed are strongly bound, they can only 
be excited above these energies. Measuring the spectra 
at these core-electron excitation thresholds to unoccupied 
states above the Fermi level bears resemblance to soft x- 
ray absorption spectra of those core electrons. This scat- 
tering process is called x-ray Raman scattering, and has 
been used successfully to measure soft x-ray absorption 
spectra of low-lying core electrons edge even in high pres- 
sure environments, accessible by the high-energy x-rays 
used in the experiment] 34 ' 35 

The IXS spectra in this regime are similar to those 
that could be studied by EELS except for the impor- 
tant possibility to access large momentum transfers (up 
to 25 nm^ 1 in this study but in principle unlimited). For 
example, an early EELS study of YBCO has revealed 
spectra equivalent to those reported here at the lowest 
q (Ref. |l2|) . The anisotropy between the ab plane and 
the c-axis results is surprisingly small in the experiment. 
Meanwhile, the theory (both within TDLDA and RPA) 
would predict a slightly larger anisotropy^ at least be- 
tween the calculated a and c directions. 

The spectral features found here are labeled with let- 
ters A-G in Fig. [T] The features A and B correspond 
to the excitation of the Ba 5p electrons (the Ba (?2,3 ab- 
sorption edges). The features C-E are rather intriguing. 
The relevant S(q, to) range for the case of q | c is showed 
zoomed-up in the Fig. [2J in both TDLDA and RPA to- 
gether with the experiment. It is in this region where 
the differences of TDLDA and RPA are the most pro- 
nounced. In RPA the peaks D and E are broad and weak, 
but appear better defined in TDLDA. On the other hand, 
experimentally there are only two excitations, which we 
label tentatively as C and D + E since it seems that the 
peaks D + E have merged together within the experimen- 
tal accuracy. All these peaks appear at higher energies 
in the experiment than in the theory. The energies of 
the excitations C—E disperse with momentum transfer, 
for the excitations D + E most strongly in the experi- 
mental result. As a guide for the eyes, their dispersion 
is marked with dashed lines in Fig. [I] In the experiment 
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FIG. 3: The real part (left) and imaginary part (right) of the 
dielectric function as calculated within the TDLDA for the 
q || c direction. The excitation dispersions are labeled as in 
Fig. [T] and nominal core-excitation thresholds are annotated 
on top. 



C and D + E seem to coalesce into a single peak at low 
momentum transfer. To observe these excitations sepa- 
rately, an access to large momentum transfers is essential, 
in particular since it seems that the peak E starts gaining 
considerable spectral weight only for q > 20 nm" 1 . These 
peaks are reproduced by the theory in both approxima- 
tions but a clear discrepancy between experiment and 
theory remains. In particular in the experiment the ex- 
citation D + E seems to follow a quadratic dispersion 
typical for a plasmon, a trend not reproduced by theory 
even if TDLDA reproduces the excitations better. This 
could be one of the key differences where a treatment of 
correlations beyond TDLDA is necessary ! 37 i 38 These ex- 
citations are sensitive to correlation, as suggested by the 
evident differences between the two exchange-correlation 
approximations. It would seem natural to assign them 
to Y 4p which has an excitation threshold in the same 
energy. However, surprisingly, by occupation analysis we 
have found them to originate mostly from Ba 5p instead. 
As the wave functions of these excited states overlap with 
those of Y 4p, the interaction of Ba 5p and Y Ap excita- 
tions seems to be strong. Finally the features F and G 
are related to the Ba 5s and Y 4s + O 2s excitations, 
respectively. 

The dielectric function as calculated within the 
TDLDA is shown in Fig. [3] for the q 1 1 c direction. Any 
single-particle excitation, associated to a transition of 
a single electron from an occupied to an empty state, 
should be apparent as a peak in Im(£r). The position 
of such an excitation corresponds to the transition en- 
ergy between the two states of the band-like picture, 
eventually renormalized by electron-electron interaction 
effects (single-quasiparticle 3 ^ excitations) and electron- 
hole (excitoni o 40 ! 41 ) effects. 

On the other hand, collective plasmon excitations 
would manifest themselves as zero-crossings of the real 
part of the dielectric function Re(e), when the imagi- 
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nary part Im(e) is small. These frequencies are in 
correspondence with the self-sustained resonant modes 
of the system, associated to collective oscillations of the 
electron plasma. When the plasmon is strongly damped 
(usually due to an interaction with single-particle exci- 
tations), the zero-crossing disappears but a relict may 
remain as a minimum of Re(e). Here we find no zero- 
crossings of Rc(e) in agreement with earlier studies^ 
implying that the plasmons in this system are heavily 
damped. While the peaks A-C are produced by dips in 
the real part at low q, they change their character into 
a peak in the imaginary part at highest q studied here, 
showing a transition from a highly-damped plasmon to 
a single-particle excitation as q is increased. The peak 
D has the clearest damped-plasmon-type character with 
the largest minimum in Re(e) at low q. While much of 
the spectral weight is due to the Cu 3d and O 2p, most 
of the fine structure is in fact due to the Y Ap and Ba bp 
electrons. By occupation analysis we have identified the 
peaks A-D to originate from Ba bp, but especially at the 
largest momentum transfers the Y Ap single-particle core 
excitations overlap with the excitation D. In other words 
the Y Ap and Ba bp core excitations interact strongly with 
the Cu 3d + O 2p valence excitations. For the peak E 
there is a weak peak in Im(e) at the corresponding en- 
ergy and momentum, suggesting a single-electron type 
excitation for this peak, mostly Y Ap character. 

A surprising aspect in the charge response of YBCO 
in this energy range is that it seems to fall outside 
the above-described normal separation between collective 
and single-particle excitations, and rather has a peculiar 
mixed character. For instance, we found that at low q the 
changes in Rc(e) due to Ba bp, often considered as a core 
state, seem to be indispensable for obtaining a quantita- 
tive agreement between the theory and the experiment. 
The Ba bp excitations introduce new minima in Rc(e) 
around 25 eV reducing its value by approximately 50%, 
giving a large contribution to the collective C-D exci- 
tations at low q. Furthermore, below 20 eV the Ba bp 
electrons increase the value of Re(e) and decrease the 
spectral weight in this energy range in the final loss func- 
tion. This is even true below 10 eV where the increase 
in Re(e) is approximately 10%. However, at q > 20nm _1 
their response is seen as peaks in Im(e), indicating in 
turn typical core-electron excitations. On the contrary, 
the Y Ap excitations have mostly a single-particle char- 
acter at all values of q studied here. Nevertheless, they 
have a significant role in determining the final shape of 
the loss function. 

The results shown here are a clear proof that the va- 
lence dynamics of YBCO cannot be described in detail 
by considering only the Cu 3d and O 2p, but at least the 
Ba bp and Y Ap have to be taken into account. Besides 
hybridizing, they couple to the Cu 3d + O 2p through 
the local field effect and they may be important also at 
much lower energies than the ones investigated here. 

All real electronic systems, and even the jellium model, 
exhibit damping of the collective excitations at increased 



wave vectors — a result well known from Landau's the- 
ory (see for instance the discussion in Refs . l42l and l43l and 
references therein). One of the reasons for the plasmon 
damping is the coupling to single-particle electron-hole 
pair excitations. In simple metals the effect of single- 
particle excitations on the plasmon lincwidth has been 
studied using, e.g., nearly- free-electron approach^ as 
well as more recently using an ab initio method^ In all 
cases where single-particle excitation channels arc strong 
in the range of plasmon energies, the damping will be 
more severe than in the jellium model. This is exactly 
the reason why a single-band theory cannot describe 
this situation: all electronic excitation pathways, which 
depend very much on the full band structure, have to 
be taken into account equally. In the cuprates and other 
so-called strongly correlated systems this may not be a 
trivial task. Even excitations of semicore electrons can 
play an important role in the damping of collective exci- 
tations, as we have shown in YBCO. On the other hand, 
it is not yet clear whether the band-like paradigm still 
holds for the highest occupied valence levels in these sys- 
tems. This was exactly the motivation for the current 
work. We have shown here that ab initio TDDFT can 
describe a complicated system like YBCO and increase 
our understanding of the interplay between semicore and 
valence electrons. 



V. CONCLUSIONS 

In conclusion, wc have presented a combined experi- 
mental and theoretical study of the dielectric function of 
YBCO for large wave vectors, probing small interparticle 
distances. This previously unexplored area of studies re- 
veals new interactions between Cu 3d and O 2p valence 
excitations with the weakly bound Ba bp and Y Ap elec- 
trons in the dynamic response and screening. Especially 
the role of Ba bp is found to be of mixed single-particle 
and damped-collectivc type giving rise to a novel exci- 
tation, observable only at large wave vectors, suggesting 
its importance in screening at short distances. We have 
shown that the valence dynamics interacts strongly with 
the low-lying semicore states. We have answered as well 
the question of how well we are able to model a HTSC 
cupratc system with many-particle interactions of the va- 
lence electrons that invalidate the band structure picture. 
In fact, against all expectations, ab initio TDDFT proves 
to be a valid and powerful theory to describe a HTSC 
cuprate and gives results which agree with the experi- 
ment very well. Finally, it is possibly in the remaining 
differences between the experiment and theory presented 
in this work where the effects of strong correlations, and 
perhaps even the origin of the Cooper pairing, could be 
looked for. 
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